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ATTITUDE DETERMINATION OF SPIN STABILIZED 
SPACECRAFT USING A DIGITAL ASPECT SENSOR 

ENERGY AND MOMENTUM CONSIDERATIONS 

The problem of the force free motion of a rigid body has been extensively 
treated by a number of people and a great many methods for solutions have been 
formulated. The most common and most useful for our discussion was formulated 
by Euler. The Euler angles are shown in Figure 1. They are defined as three 
successive rotations performed in specific sequence. The sequence is as follows: 
(1) Rotate the initial system of axes counter-clockwise about the Z axis by an 
angle 0, which produces a set of intermediate axes (x', y', z' ). (2) Next rotate 
the intermediate set of axes (x', y', z' ) counterclockwise through an angle 6, 
about the x' axis forming the line of nodes. This produces a second set of inter- 
mediate axes (x", y", z"). (3) Finally, rotate this set of axes counterclockwise 
through an angle 0 about the z" axis, forming the (x, y, z) axes. 

The three angles 0 , </>, and 0 then completely describe the orientation of 
the x, y, z coordinate system with respect to the X, Y, Z coordinate system. 

Now let the X, Y, Z coordinate system define an inertial space. Assume 
a rotating satellite situated in this inertial space where the momentum vector, 

(L), of the rotating satellite is orientated along the Z axis of the inertial space, 
x, y, and z may now represent a coordinate system fixed in the rotating sat- 
ellite. The axes x, y, z are aligned with the three principle moments of inertia 
of the satellites I r I 2 , and I 3 . The instantaneous values of momentum about 
the x, y, and z axes become p x , p y , and p z respectively. 

P = 1, w 

1 X 

p = I, &) 

y 2 y 


p z = 1 3 W Z 


( 1 ) 

( 2 ) 

( 3 ) 


1 



r 


z 



Where o; y , and co x are the instan- 
taneous values of angular velocity about 
x, y, and z. 

From examination of Figure 1 it can 
be seen that the components of L are: 

P x = L sin 0 sin i/; (4) 

P y = L sin 9 cos \p (5) 


Figure 1-Euler rotation angles 


p = L cos 9 

r z 


( 6 ) 


By the process of infinitesimal rotations 1 a relationship is established 
between the angular velocities co x , co y , and co z , and the rate of change of the 
Euler angles ft, ft, and 9. These relations are stated in equations (7), (8), and 


O). 


oo* - ft sin 6 sin ft + 6 cos ft 

co^ = ft sin 6 cos \p - 9 sin \p 

oo* = ft cos 6 + ft 

Substituting 7,8, and 9 in 1, 2 and 3 gives 

P x = Ij (eft sin 6 sin ft + 9 cos ft) 

p y -- I 2 (ft sin 6 cos ft - 6 sin ft) 

p = I 3 (ft cos 9 + ft) 


(?) 

( 8 ) 

( 9 ) 


( 10 ) 

( 11 ) 

( 12 ) 


^Pg. 134, Goldstein, H. - Classical Mechanics, Addison-Wesley Publishing Co., Reading, Mass. 
1957. 
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Multiplying (4) by sin i p and (5) by cos \p gives, 

P x sin 0 = L sin 6 sin 2 ip (13) 

P y cos ip - L sin 6 cos 2 i p (14) 


Now add expressions (13) and (14) to get 

L sin 6 - p x sin i// + p cos \ p (15) 

Assume for simplicity that the satellite is balanced in such a way as to make 

= V 

Making this assumption, and substituting (10), and (11) in (15) gives 

L sin Q - Ij <f> sin 6 (16) 

or 

L = Ij 4> (17) 

Again, multiplying (4) by cos 0 and (5) by sin 0 and subtracting gives the 
following expressions. 

P x cos ip - p y sin ip - L sin 6 sin ip cos ip - L sin 6 sin ip cos ip (18) 

p x cos ip - p y sin ip - 0 (19) 

Substituting expressions (10), and (11) into expression (19) gives 

0 = 0 ( 20 ) 

or 6 is time independent. 
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Substituting expression (17), into expression (6), and then (6) into (12), 


Ij 0 COS 0 = I (<j> COS 8 + 41 ) 


( 21 ) 


or, 


a) 



0 cos 8 


b) 


4 > = 



(Ij - I 3 ) COS 6 


( 22 ) 


\jj is the angular velocity of the satellite about the satellite z-axis. This 
will be called spin rate. This is to be distinguished from the apparent rotation 
rate of the satellite with respect to a fixed external point. This apparent rotation 
rate has an average value of (0 + 0 ). 

0 is the angular velocity of the satellite z-axis about the momentum vector 
(or more precisely, the angular velocity of the line of nodes). This will be called 
the precession rate.* 

Q is the rate of change of the precession half-cone angle. Such a motion is 
called nutation. * A satellite with appreciable nutation presents an extremely 
difficult, if not insoluable, motion determination problem. Nutation is the result 
of I j / 1 2 and 0 / 0 . 

As previously mentioned, if Ij is made equal to I 2 , 8 = 0 and there is no 
nutation. 

0 is the precession cone half angle. 

A freely rotating satellite which has both spin and precession motions ex- 
periences time dependent accelerations at various points in its structure. Since 
most satellites are not absolutely rigid bodies, some flexing of the structure 
occurs and energy is thus dissipated through frictional losses. Momentum is, 


*NOTE: The use of the terms precession and nutation are defined in a force free field according 
to the astronomers’ convention, i .e. : Precession rate is the rate of rotation of the lines of nodes 
and nutation is the short term fluctuation in the precession cone half angle. This convention is 
quite different from that used in the theory of gyroscopic motion. Gyroscopic precession is the 
motion of the gyroscopic spin axis as a result of an external torque applied perpendicular to the 
spin axis. Gyroscopic nutation is the wobble of the gyro axis about the momentum vector. Nuta- 
tion in the gyro sense is similar to precession in the astronomical sense. 


4 



of course, conserved as long as no external forces act on the satellite. The 
effect of such energy losses can be seen by examining the following equations: 

The Lagrangian is written 


L = T - V = 1/2 (I, col 


+ I 2 col + *3 > - V (d ' ^ 


(23) 


Since we are considering force free motion, V (0, 0, 0) = 0. 

Substituting (7), (8), (9), in (23) gives 

T = 1/2 [ij (0 sin 0 sin 0 + 0 cos 0) 2 + I 2 (0 sin 6 cos 0-0 sin 0) 2 


+ I 3 (0 cos 0 + 0) 2 ] 


(24) 


Assuming again that the satellite is balanced so that I t = I 2 , expression 
(24) reduces to 


T = 1/2 I, (0 2 sin 2 0 + 0 2 ) + 1/2 I 3 (0 cos 0 + 0) 2 (25) 


Since 0 is a constant (equation 19) , 0 = 0. Now substituting equation (22a) 
in (25) gives 


1/2 Ij 0 2 ^sin 2 0 + — i cos 2 0^ 


(26) 


Substituting L = Ij 0 from (17) gives 


L 2 

21 , 


T - (sin 2 0 + — cos 2 8^ 


(27) 
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In the absence of any external forces L is a constant. If precession motion 
causes energy to be dissipated internally through battery liquids sloshing, non- 
rigid members flexing, etc., T becomes smaller. Since L is not affected by 
internal dissipation of energy, the value of Q then must change; as shown in Fig- 
ure 2, the value of Q will tend toward 0°, or 90°, or remain constant depending 

Case I - If I j /I 3 < 1, 9 tends toward 
0°. In this case the precession-coning 
will damp out in time. Since it is usually 
desirable for the satellite to rotate about 
the z - axis, most spin stabilized satellites 
are balanced so that the z -axis coincides 
with the largest inertia moment, that is so 
that Ij/Ij < 1 . 

It should not be assumed that damping 
of the precession cone implies that <p 
goes to zero. Equation 22 states that when 
6 = 0 °, 4> = I 3 4>/I 3 - I 3 . It can be seen 
from Figure 1 that the apparent rotation 
velocity of the satellite for 6 = 0 ° is 
+ \p )• 

Case II - On the other hand if 
Ij/I 3 > 1 , 6 tends toward 90°. In this 
case the precession cone opens up until 
the satellite begins to spin like a propeller. 
In this mode 0 goes to zero (Equation 22 ). 

Case III — The third possibility is that 
I /I 3 = 1. In this case the satellite is 
perfectly symmetrical about all axes and 
there is no way of distinguishing the z - 
axis dynamically. In this case I = I 3 
and from equation 22 again <jj = 0 . 

The final result of any of these three cases is to cause the satellite to rotate 
with no wobbling motion about the axis of the largest moment of inertia. 

It should be noted that Case III is difficult to achieve in reality. 1 x is sel- 
dom made exactly equal to I 3 and in fact, our symmetry assumption that I t = I 2 
is seldom precisely true either. Thus, it may happen that if the three moments 


on the value of the ratio I /I 3 . 



0 CASE I < I 8 TENDS TOWARDS 0° 




Figure 2-Relationship of energy to preces- 
sion half angle for constant momentum. 
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are very nearly equal, the condition where Ij > I 3 > I 2 , orI <I 3 <I 2 may 
inadvertently occur. This possibility raises a fundamental problem associated 
with asymmetrical rotating bodies.* That is, if a body is shaped so that either 
Ij > I 3 > I 2 , or Ij < I 3 < I 2 , and the body is set rotating about I 3 , its motion 
is unstable. This instability is such that any slight perturbation will send the 
body into wild, unpredictable, and quite violent gyrations. The change in spin 
axis caused by this type of instability occurs in a time interval on the order of 
one revolution period and requires no energy dissipation. Thus, this instability 
is quite different from Case II where the change in spin axis occurs over many 
spin periods and depends entirely on the dissipation of energy. 


MEASUREMENTS USING THE DIGITAL ASPECT SENSOR 

To make the measurements necessary to determine the motion of a spinning 
satellite, assume that a photo detector is mounted behind a slit situated on the 
side of the rotating satellite. Let the axis of this detector be along the satellite 
x -axis with the slit positioned parallel to the satellite z-axis. Then the field 
of view of this detector will be a fan which will sweep across the sky, which is 
the celestial sphere, as the satellite rotates. This shown in Figure 3. Now 
further assume that this fan field of view is quantized into a number of segments 
as shown in Figure 4. This may be done if instead of a simple slit and photo 



SATELLITE 
2 AXIS 



CELESTIAL 
SPHERE 


Figure 3-Celestial sphere representa- 
tion of a precessing satellite using a 
digital aspect sensor. 


Figure 4-Quantized field of view of a 
digital aspect sensor 


NOTE: For a rigorous treatment of the stability of asymmetrical rotating rigid bodies see page 
116 of "Mechanics” by Landau, L. D. and Litshitz, E.M., Addison-Wes ley Publishing Co., I960. 
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detector, a digital aspect sensor 3, 3 is used to provide the fan shaped field of 
view. Each of the fan segments has a unique identifying number. As the satellite 
rotates, if a bright object such as the sun should be in a portion of the sky swept 
out by one of these quantized fan segments, the number corresponding to that 
particular segment would be transmitted to the ground. This number then defines 
the angle between the satellite z -axis and the sun vector. This angle must lie 
between the upper and lower edges of the particular quantized fan segment. The 
angles between the edges of these quantized segments and the z -axis are listed 

in a calibration table and are designated a y , a 2 , a 3 a n . If 

there is no precession coning, the same segment number will be repeatedly 
transmitted upon each revolution of the satellite. 

Most of the discussion in the rest of this section deals with problems caused 
by precession coning. For the zero precession case the spin rate 0 and preces- 
sion rate 4> can be found from equations (38) and (39). The reader may then skip 
to equations (51) and (52) where the initial phase angles \p and <p 0 can be set 
equal to zero at time t 0 . Much of the following discussion of precession motion 
determination will reduce to a degenerate case for zero precession. 


Calculation of fi and 0 


In the event of precession coning the motion of the fan across the celestial 
sphere is quite complicated. The first quantities to be calculated from the data 
are the angles fi and 0 . 

If the satellite is balanced so that Ij = I 2 the angle a between the z-axis 
of the satellite and the sun vector will be a function of time given by (28) 

cos a - cos fi cos 6 + sin 0 sin fi sin <p t (28) 

Where P is the angle between the momentum vector and the sun vector. Figure 
5 is a graph of (28) showing sun sightings for the two cases where I t > I 3 , and 
Ij < I 3 . If I > I 3 , the fan "sees” the sun many times for one cycle of <f>. If 
Ij < I 3 , 4> completes more than one cycle per sun sighting. Often the best 


3 Albus, James S., "A Digital Solar Aspect Sensor,” Astronautics, January 1962 also NASA TN 

D-1062. 

O 

°Albus, James S. and Schaefer, David H., “Satellite Attitude Determination: Digital Sensing and 
On-Board Processing,” IEEE Transactions on Space Electronics and Telemetry, Vol. SET-9, pp. 
71-77, September 1963. 
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Figure 5-Typical plots of sun sightings from a 
digital aspect sensor 


method of calculating /3 and 6 is to plot the a limits for the sun sightings and 
fit a curve. It is possible to apply a statistical approach to the calculation of fi 
and 6 over a period which is long compared with any beat frequency set up be- 
tween 4> and 0 . However, this is sometimes in practice a prohibitively long 
period of time. 

To calculate the values of /3 and 6, take some time interval Tj - T 0 such that 
at Tj and T Q (0 - v/2) < and (V' + tt/2) < where and g 2 are arbitrarily 
small angles. From Figure 3 it can be seen that these restrictions make a a 
minimum at T and T . Between T 0 and Tj the sun appears at semi-regular 




intervals which are not integrally related withc^. Therefore, over a period which 
is long compared to any beat frequency which might appear between and 0, all 
values of (*j>t) between 0 and 2tt become equally probable. 

Now the probability of a sun appearance in the segment nearest the spin axis 
is equal to N a /N (where N a is the number of appearances of the sun in the nearest 
segment, and N is the total number of times the sensor fan crosses the sun during 
Tj - T Q ). Thus, the probability at any instant of a > a a is N a /N. 

A half cycle of the curve (28) is plotted in Figure 6. 


a 



Figure 6-Plot of 1/2 cycle of 
the curve defined by Equation 28. 


If an arbitrary point along the 4> t axis is 
chosen, the probability that the corresponding 
value of a is less than a is 2t /T,. The re- 

a a <P 

lationship between a a and t a is given by sub- 
stituting in (28) 


cos a g = cos ft cos 9 + sin / 3 


sin 6 cos (<£> t a ) 


(29a) 


which can be expressed as 


/ 2 t a \ 

cos a g = cos ft cos 6 + sin ft sin 0 cos I 77 — — 1 


(29b) 


We have already shown that the probability of a solar appearance during the 
interval Tj - T Q in the nearest fan segment is N /N. If a g is set equal to the 
value of the most distant edge of the nearest quantized segment, then N a /N may 
be substituted for 2t a /T in (29b) to give (30). 


cos e cos /3 cos 9 + sin ft sin 6 cos 



(30) 


Similarly, 


cos a b = cos ft cos 6 + sin ft sin 6 cos 



(31) 
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where a g is the most distant edge of the nearest quantized fan segment to the 
spin axis, andN a is the number of sun appearances in the nearest segment. a b 
is the nearest edge of the most distant quantized fan segment. N b is the number 
of sun appearances in all segments except the most distant. N is the total number 
of all sun appearances during the interval Tj - T 0 . 

The restrictions on T Q and Tj may now be disregarded by making N sufficiently 
large. 

Writing (30) and (31) in more convenient form for solution gives 


cos (/ 3 + 6) 


ad-bc+a-c_a(d + l)-c(b + l) 


d - b 


d - b 


cos (yS - 6) - 


ad - be - a + c a (d - 1) - c (b - 1) 


d - b 


d - b 


(32) 


where 


a = cos a 

a 

b = cos (?7N a /N) 
c - cos a L 

b 

d = cos ( 7 TN b /N) 

Note that equation (32) is ambiguous as to the quadrant of (/3 + 6) and {fi - 6). 

Thus the digital sun sensor does not indicate directly whether the sun is inside 
or outside of the precession cone. Ignorance of this ambiguity can lead to con- 
siderable confusion if the sun happens to lie inside the cone. Resolution of the 
ambiguity is not straightforward unless the approximate size of the precession 
cone can be independently measured. 

If the rotation rate can be determined by an independent means such as a 
magnetometer or a telemetry AGC record use of equations (36) and (37) will 
provide an indication of whether the sun is inside or outside the precession cone. 
However, if the magnetic field vector or the ground receiving station is inside 
the precession cone this method itself is subject to the same type of ambiguity. 
Another method is to reconstruct a time history of a precession cycle on the 
celestial sphere. The phase of the 4> and \p angles will usually provide an indi- 
cation of whether (/3 - 6) is plus or minus. 
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Where the sun only appears in two quantitized segments, there are not 
enough known quantities to solve for both /3 and 6. In this case, unless 9 can 
be estimated from previous data, assume 6 equal to width of smallest segment 
in which the sun appears. Then solve (33) for / 3 . 


cos (i - 


ab + c / (b 2 + c 2 ) - a 2 
(b 2 + c 2 ) 


(33) 


0 < /3 < 180° 


where 

a = cos a 

a 

b = cos 6 

c = sin 6 cos (tt N 0 /N) 

In this case, if 6 is constant over a period of days, an accurate measure of 
6 can sometimes be made by noticing how long the sun's appearances oscillate 
from one segment to the other. The velocity of the sun's yearly motion, multi- 
plied by the crossover oscillation time, give 2 9. /3 at the midpoint of the cross- 
over time is exactly the value of the calibrated angle a n . Other factors besides 
the sun's apparent motion can of course cause to slowly vary over a period of 
days. These factors are external torques caused by magnetic fields, radiation 
pressure, aerodynamic drag, plasma winds, gravitational gradients, etc. 


Calculation of <p and 0 


Having now calculated 6 and /3 we proceed to determine the spin rate i/ 1 * 
and the precession rate 4 > . In the general case where precession motion is 
present, the apparent rotation period, which is the time between successive 
crossings of the sun by the fan field of view, will not be uniform. This is due to 
the fact that sometimes the fan is tilted forward when it crosses the sun, and 
sometimes backward. However, the average value of the apparent rotation period 
is constant. In fact, and i p may be determined from the average apparent 
rotation period by the following method: 
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Choose X in Figure 3 so that the sun lies in the XOZ plane. This may be 
done with no loss in generality. Choose a time that the fan crosses the sun such 
that (cf> - 77/2 ) < and (i p + v/2) < g 2 where ^ and are arbitrarily 
small angles. Call this time T 0 . This condition can always be satisfied at some 
time since 0 and 0 are not integral multiples of each other. Now wait for another 
time when the fan crosses the sun such that these same restrictions can be ap- 
plied to 0 and 0 . Call this time T x . Count the number of times the fan crossed 
the sun between T Q and T x . Call this N f . 

Then the average period 


T 

a v 


T 


0 


N 

r 


Now from T Q til T x , 0 has made 
number of degrees). Therefore, 


revolutions (within an arbitrarily small 


2^ = 0 (T x -T 0 ); 


(34) 


and during that same time 0 made revolutions. Therefore, 


2^=0CT x -T 0 ). 


(35) 


Now if Ij > 1 3 the value of can be immediately determined by inspection 
of a plot similar to Figure 5. In fact in the case where I 3 > I x , if one has a 
rough knowledge of Ij/I 3 > can usually be found by inspection. However, if 
there is any doubt about , or if the sun lies inside of the precession cone a 
more rigorous approach is more reliable. 

Assume for our discussion a fan field of view 180° in length and negligibly 
thin. In Figure 4 this would mean that a x = 0 and a n = 180°. From Figure 3 
we can write 


N r =% + N * i f 


and 

N r > if p<e 


(36) 


(37) 
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The case where N r > in equation (37) arises when the sun lies inside 
certain regions within the precession cone. These regions are quite complicated 
and will not be discussed here. 

For the case where the fan is not 180° in length (i.e. a x > 0 and/or a n < 
180°) equation (36) holds where (6 , + a 1 )< /3 < (a n - 9) and equation (37) holds 
where (0-a.j) > J3 > [360° - (a n + 0)]. If these conditions are not satisfied, a 

value for N f can sometimes be determined by inspection of a plot similar to 
Figure 5. 

If the conditions for equation (36) hold, it can be combined with equations 
(22), (34), and (35) to give equations (38) and (39). 


<t> = 

(T, - T 0 ) 


2ttN 

r 



(38) 


and 


0 = 


(T, - T 0 ) 


2 7tN 

r 



cos 6 


(39) 


Now if the restrictions on choosing the times T 0 and T are discarded, the 
uncertainty in (0 + 0) will be less than ± (50%/N r ). Therefore, if T 0 and T 
are merely chosen a large number of revolutions apart so that N > 50 (38) 
and (39) give quite accurate answers. 

The accuracy of 0 and 0 from equations (38) and (39) depends, of course, 
on the accuracy of Ij/I 3 and 9. If 0 and 0 can be roughly determined from 
equations (38) and (39), and then the knowledge of 0 and i p used to find by 
inspection, equations (38) and (39) can then be solved again with 0 and 0 as 
known and Ij /I 3 unknown. This kind of iteration can be used to refine the value 
of I j /I 3 and even under certain conditions to determine 6 . 


Calculation of 0 Q and 0 o 

The integrals of motion 0,0, 9 having now been calculated, it is necessary 
to establish initial conditions in order to define the instantaneous position of the 
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satellite in inertial space at a particular time. A particular passage of the slit 
across the sun may now be chosen as occurring at t 0 . 

In Figure 7 the digital fan is shown crossing the sun. a a and a fa correspond 
to the angles between the spin axis and the upper and lower edges of the quan- 
tized fan segment containing the sun. From this figure the following equations 
may be written 

cos a = cos (3 cos 6 + sin (3 sin 6 cos (90° - 4>) (40) 

or in the case of a and a, 

a b 

cos a a = cos fi cos 6 + sin /3 sin 9 cos (90° - <£ A ) (41) 


cos a b = cos /3 cos 0 + sin/3sin0cos (90° - <P B ) (42) 

Substitute the relation 

4> - 0 (t - t 0 ) + 4> 0 . (43) 

cos a a = cos /3 cos 6 + sin /3 sin 6 sin (4> ( t - t Q ) + 0 OA ) (44) 
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At t = t 0 (44) and (45) can be solved for <£> 0A and 4> QB . 4> 0 must lie some- 
where between <p QA and 4> 0B . However, the solution of (44) and (45) is ambiguous 
regarding the quadrant of 0 OA and 4> QB . This ambiguity may be resolved by 
choosing a value of t at a sun sighting approximately an odd number of quarter 
periods of 4> away from t Q . Solve (44) and (45) again for cos a a , and cos a fc . 

This time try first one possible quadrant value for 0 QA and <t> 0B and then the 
other. Compare the two values of cos a a and cos a fa obtained with the cosine 
of the actual value of a and a L measured at time t. The closest value should 
indicate the proper quadrant for <£ QA and cf> . 

Having now determined unique values for <fi QA and cp QB , it is possible to 
narrow the limits on the value of 4> 0 . The precession rate 4> is accurately known, 
and the time of each sun sighting is also well known. Therefore, equations (44) 
and (45) can be repeatedly solved for 0 OA and 4> QB by inserting <x & , a b , and t 
for each sun sighting after t Q . Since 0 O must always lie between 4> OA and <fi OB , 
4> 0 must lie between the maximum of the smaller values and the minimum of 
the larger values of the pair <p QA and 4> OB . This is somewhat tricky when put 
into practice since whether <fi QA or 4> QB is the larger of the pair depends on the 
quadrant of the argument [<j> (t - t 0 ) + 0 O ] . If cos [<f> ( t - t 0 ) + <£ 0 ] is positive, 
^oa > ^ob » but if cos [0(t - t 0 ) + 0 O ] is negative, <£ QA < 4> OB . 

Once 4> 0 is determined accurately enough, solve (47) for cos a Q . 

cos a Q = cos J3 cos 9 + sin /3 sin 6 cos (90° - <fi Q ) (47) 

where 0 < a Q < 180° (48) 


Now from the same spherical triangle in Figure 3 


COS 0 O - 


sin /3 cos 0 O 


(49) 


sm a„ 


sin 0 O = 


cos /3 - cos 9 cos 
sin 6 sin a. 




(50) 


Solve (49) and (50) for 0 O . 
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The instantaneous value of t p and 0 can now be determined by 


0 = i> (t - t 0 ) + 0 O (5i) 

4> = 4> (t - t 0 ) + 0 O (52) 


These equations will hold for some length of time t - t Q . When errors 
begin to accumulate, a new t Q should be chosen and a new set of initial conditions 
computed. 

If the time base in the spacecraft is quantized, there is a quantizing under- 
tainty in the value of t Q . Equations for removing quantizing uncertainties in 
the value of t 0 under the conditions of 6 = 0 are contained in Appendix E. 

MEASUREMENTS USING A HORIZON DETECTOR 

Thus far all the measurements have been made on the sun. Since the sun is 
considered a point source at infinity, there is an axis of symmetry from the center 
of gravity of the satellite to the sun. Another inertial reference point is necessary 
in order to eliminate the ambiguity caused by this symmetry. In general what 
is needed is another known vector which is not colinear with the sun vector, and 
an angle § which is measured from this second vector to the spin axis. This 
second vector may correspond to the magnetic field, to the subsatellite zenith 
vector, or any other desirable vector. 

In the following discussion the use of horizon sensors to measure an angle S 
between the spin axis and the subsatellite zenith is considered. However, the 
equations in Section IV apply as well to a 8 angle measured by any other type 
of sensor from any other type of second vector. 

Infra-red horizon sensors which detect the temperature discontinuity caused 
by the warm earth against the cold background of outer space can provide horizon 
sensor data on the dark side of the earth as well as on the sunlit side. Visible 
light horizon detectors can be used when data is needed only on the sunlit side. 
Conditions when visible light detectors can be used are in Appendix F. Either 
type detector usually has a pencil beam field of view arranged so as to provide 
a pulse upon crossing either into or out of the earth's disc. Sometimes it is 
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necessary to transmit the knowledge of whether a pulse represents an "in" or 
"out" crossing, but for satellite applications the difference can usually be deter- 
mined by relative pulse spacings or prior knowledge about the approximate spin 
axis orientation. 



Figure 8— Definition of variables 
needed to solve the horizon sensor 
problem 


Information is contained in the repe- 
tition rate of the pulse pairs, in their 
spacing, and in their position relative to 
the sun appearances in the quantized fan. 

The repetition rate of the pulse pairs can 
provide a measure of the rotation rate. 

The information contained in the pulse 
spacing provides a measure of 8 , which 
in this discussion is the angle between the 
spin axis and the subsatellite zenith 
vector. The subsatellite zenith vector is 
a vector from the center of the earth through 
the center of gravity of the satellite. 

If the earth is considered to be a disc 
inside the celestial sphere, where the 
center of gravity of the satellite is the 
center of the sphere, an expression can be 
derived relating the pulse spacings to 8 . 


Figure 8 represents a satellite rotating with no precession so that 8 - 0 . 

An earth horizon pipper is mounted at an angle from the satellite Z - axis so as 
to sweep out a cone which cuts the earth horizons at A and C. From Figure 8b 
it can be seen that 


EB 

EA 



(53) 


where At A is the time between horizon pulses and T is the rotation period. 
From the figure, various other relationships can be derived: 


EB = ED + DB (54) 

ED = OE tan (90° - 8 ) where OE=OHcosy (55) 
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DB = OG cos (90° - S) where OG = OH cos p 


(56) 


EA = OH sin 7 


(57) 


Substituting all these values in expression (53), gives 


OH cos p • 

OH cos 7 cot S + = OH sm 7 cos 

sin S 



which reduces to 


cos p + cos 7 cos 8 = sin y sin 8 cos 



where 


(58) 


(59) 


sin p 


r E + h 


(60) 


where 

r_ = mean radius of earth 

E 

h = apparent altitude of satellite 
Rewriting so as to solve explicitly for S 

t -de ± f '/(e 2 + f 2 ) - d 2 

COS 0 

(e 2 + f 2 ) 

0 < S < 180° 


( 61 ) 
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where 


d = cos p 
e = cos y 

f = sin y cos (wAt A /T) 

Expression (61) is plotted in Appendix A. It double valued for most of the 
important cases. In order to minimize this ambiguity it is usually best to make 
y = p or 180° - p. In order to remove the ambiguity altogether a second 
horizon sensor mounted with a different y angle can be used. If two sensors 
are used it is usually best to mount one at an angle y g = p and the second at 
an angle y b = 180° - p. With a system using two such sensors, the absence 
of a signal from one or the other of the sensors usually will eliminate the am- 
biguous solution from (61). When both sensors scan the horizon the solution for 
8 becomes single valued, and in addition the satellite apparent altitude can be 
calculated. 


77 At 0 

cos p + cos y g cos 8 r sin y g sin 8 cos — — — (62a) 

77 At. 

cos p + cos y b cos 8 - sin y b sin 8 cos (62b) 


where At a is the pulse spacing from sensor A and At is the pulse spacing 
from sensor B. Subtracting (62b) from (62a) and reducing gives 


tan 8 = (63) 

1 - m 

0 < S < 180° 


where 

p = cos y a 
q “ cos y b 

1 = sin 7 a cos (77"At /T) 
m = sin '/ b cos (77At b /T) 
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Solution of equation (63) is plotted in Appendix B. 

Substituting the values of S just calculated back into equation (62a) and 
solving for p gives 


cos p = sin y a sin 8 cos ( — - cos 7 a cos § 


( 64 ) 


° < P<^ 

2 


Rewriting equation (60) gives the apparent satellite altitude 


h = 



(65) 


This apparent altitude may differ from the true altitude by the thickness 
of the atmosphere. 

Information determining the value of § is also contained in the relative 
position of the horizon pulses to the quantized fan crossing of the sun. If the 
satellite is spinning so that 6 = 0 °, horizon pulses will be symmetrical about 
the time at which the quantized fan plane crosses the center of the earth. Fig- 
ure 10 illustrates the case where the horizon sensors are mounted so that their 
fields of view lie in the plane of the quantized fan. The center of the earth is 
located on the celestial sphere at (RA z + 180°) and (-D z ), where RA z is the 
right ascension and D z the declination of the subsatellite-point zenith vector. 
(180°-77) is the great circle arc from the sun to the earth's center on the celes- 
tial sphere, where 17 is the great circle arc from the sun to the subsatellite 
zenith vector. At c is the time midway between horizon pulses minus the time 
the quantized fan crosses the sun. The relationship between 8 and At c can now 
be written from the spherical triangles in Figure 9. 

cos (180° -7 i) - cos /3 cos (180° - 8) + sin /3 sin (180° - S) (66) 


cos l(4> + </) At c J 
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or 


cos t] = cos fi cos 8 - sin /3 sin 8 cos 


[(0 + 0) At c ] 


(67) 


Solving for 8 gives 


s rs ±u/(s 2 + u 2 ) - r 2 
cos 8 = 

(s 2 + u 2 ) 


( 68 ) 


where r = cos 77 
s = cos /S 

u = -sin /? cos ( 277 -At c /T) 
0 ° < 8 < 180° 
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and where 


cos ri - sin D sin D + cos D cos D cos (RA - RA ) 

1 z s z s x s z 


(69) 


where 

D = declination of subsatellite zenith 

z 

D = declination of sun 

S 

RA z = right ascension of subsatellite zenith 
RA = right ascension of sun 

S 


01t)< 180° 

Calculation of S by (68) should result in 2 solutions. This twofold ambiguity 
may be resolved by choosing the solution of (68) which most closely corresponds 
to values of 8 obtained from (61) or (63). 


Calculation of 8 is often in practice a difficult task owing to the nature of 
the horizon. Pulse pairs do not always indicate the true horizon. In addition 
the apparent altitude of the satellite must be measured not from the earth's sur- 
face, but from some vaguely defined part of the atmosphere which is often depend- 
ent on weather conditions. Finally the earth is not a perfect sphere. 


Calculation of 8 using equation (63) eliminates consideration of the apparent 
altitude of the satellite, and has proven in practice to usually provide the most 
reliable method of calculating. Use of equation (68) also eliminates consideration 
of apparent altitude. 


Perturbation analysis of the solution for S indicates that under certain con- 
ditions one type of solution for 8 is less sensitive to measurement errors than 
another. This type of perturbation analysis is perhaps most understandable when 
done graphically. Appendix A, B, and C contain plots of the solutions of equations 
(61), (63), and (68) under various conditions. Sensitivity of 8 to perturbations 
in various other parameters can be examined on these plots. 
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If there is some precession motion present so that 6 ^ 0 , some sort of 
averaging procedure is needed to obtain an average 8 . This average 8 then 
will be the angle between the second known vector and the satellite momentum 
vector. 

Sufficient information has now been obtained to calculate the instantaneous 
orientation of the satellite for any moment of time. Thus, knowing the mounting 
arrangements of any sensor on the satellite, it is possible to calculate the sensor 
orientation for any time. 


SPACECRAFT ATTITUDE DETERMINATION 

The spacecraft momentum vector position can now be calculated. From 
Figure 10 we may write equations for the values of the auxiliary angles Z 1 and 
Z 3. 


cos Zl 


cos [i - cos Tj cos 8 
sin 7] sin 8 


(70) 


where 


sin Zl 


sin P> sin Z2 
sin -q 


Z2 = (<p + 4>) At c ± 180° 


sin D - cos t? sin D 

cos Z 3 = 

sin 7] cos D z 


(71) 


(72) 


s in Z3 


sin (RA - RA ) cos D 

' s z s 

sin rj 


(73) 


Now using the values of Zl and Z3 from (70-73) we may find RA L and D L 
the right ascension and declination of the momentum vector. 
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sin D = sin D cos S + cos D sin S cos (Z3 + Zl) 

L Z Z 


( 74 ) 


- 4r < D. <1 


7 r 

2 < Ll ~ 2 


cos - RA z ) 


cos S - sin D sin D T 

z L 

cos D cos D t 

z L 


(75) 


sin (RA l - RA^ ) 


sin (Z 3 + Zl) sin 
cos D. 


(76) 


RA l = RA z + (RAj^ - RA z ) 


(77) 


This solution is unique. If the quadrant of (<£ + </) At c is unknown, Zl must 
be found from equation (70) and two solutions result. 


Error Calculations 


It is well to question at this point the accuracy of the position of the momen- 
tum vector determined by these methods. Tables of values of the error in the 
position of the momentum vector caused by a one degree error in 8 is contained 
in Appendix D. From these tables certain facts become obvious. 

The error is large whenever r\ ~ 0 , or 77 ~ 180°. The error is also large 
whenever 77 i /? + 8 , 8~77 + /3,/3~S + t 7 , or t] + {3 + S ^ 360°. 

Also contained in Appendix A, B, and C are a series of plots of S vs. At , 
At b , and At under various conditions. From these plots then it is possible to 
estimate the error in momentum vector position caused by an error in the meas- 
urement of At , At , or At c . From the plots find the slope of the curve under 
the conditions in question, multiply it by the uncertainty in At a , At b , At c . 

This gives the uncertainty in S . Then go to table and look up the error in mom- 
entum vector position corresponding to the error in 8 . 
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The table in Appendix D is symmetrical and therefore the error in momentum 
vector position caused by an uncertainty in /3 can also be obtained by interchang- 
ing the values of (i and 8 and assuming that S is the uncertain quantity. Thus 
the uncertainties from (i can then be combined with those from 8 to give the 
total uncertainty. 


SPACECRAFT 



Figure 11— Position of a sensor 
Sj relative to the aspect sensor 
in a spacecraft coordinate system 


Sensor Position as a Function of Time 


The momentum vector having been deter- 
mined, the constants of motion calculated, 
and the initial conditions at time t Q measured, 
it is now possible to reconstruct the motion 
history of the spacecraft both forward and 
backward in time from t Q over an interval 
as long as the uncertainties in the above quan- 
tities will permit. The look direction of any 
sensing device on the spacecraft can then be 
found at any time if the sensor position rela- 
tive to the aspect sensor is known. 

Assume a sensor S; mounted on the 
spacecraft so that its axis makes an angle 
ARC with the I 3 inertia axis, which is the 
satellite z-axis, and so that the sensor axis 
lies at angle AZSEN. from the aspect sensor 
axis measured east along the satellite equa- 
tor. See Figure 11. 


The following solution for the instantaneous position of the individual sensors 
does not use the conventional coordinate transformations by matrix multiplication. 
The matrix method is perhaps more elegant, but in general the number of com- 
puter operations for the matrix method is considerably in excess of that needed 
for the method described below. 


NORTH 



Figure 12— Relative positions of a sensor 
axis, the spacecraft z-axis, and the mo- 
mentum vector 
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NORTH 



Figure 13-Relationships of sensor axis Sj with 
respect to the spacecraft z-axis and the momentum 
vector 


Define an angle A as shown in Figure 12. 


A =Q+^-0 


(78) 


where </>=<£(t-t o )+0 o as defined in equation (52) and where 


sin Q 


sin (RA -RA t ) cos D 

v s L' s 

s in /3 


(79) 


sin D - cos sin D, 

Q s / L 

= ; 

sin /S cos D l 


(80) 


Also from Figure 12 


B i = (f - - azsen. 


( 81 ) 


where <//=</> (t - t 0 ) + \p 0 as defined in equation (51) 
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Now solve for the auxiliary angles E. and C. from Figure 13. 
If 6 = 0 

E. = ARC. and C. = 77 - B. 


However, if 6 > 0 


cos E. - cos 6 cos ARC. + sin Q sin ARC. cos B 
1 1 1 1 


0 5 E. <77 

i 


sin B. 

sin C. - sin ARC 

1 sin E. 1 

i 


cos C 

i 


cos ARC - cos 9 cos E 

i i 

sin 6 sin E. 


Finally the right ascension RA. and declination D. of the sensor S 
be found by 


sin D = cos E. sin D. + sin E. cos D, cos(A-C) 

1 1 L 1 L 1 



< 
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~2 


where D L is the declination of the momentum vector from equation (74), 


cos (RA. - RA l ) = 


cos E - sin D l sin D. 

1 

cos D r cos D. 

L l 


sin (A - C. ) sin E. 


( 82 ) 

(83) 

(84a) 

(84b) 

can 

1 

(85; 

and 

(86 


sin (RA. - RA ) 

1 L 


AC D 


(87 



RA. = (RA. - R\) + RA l 


( 88 ) 


where RA L is the right ascension of the momentum vector from equation (77) 

Thus the aspect problem is uniquely solved for any instant of time. It is 
now possible to obtain the angle between the sensor S. and an external vector 
such as the magnetic field vector or velocity vector. The coordinates of the S. 
vector need only be substituted into equation (89) along with the coordinates of the 
external vector to obtain 


cos Q. = sin D. sin D + cos D. cos D cos (RA. - RA ) 

1 1 X 1 X ' 1 X 


(89) 


0 1 fi. < n 

1 


where a = the angle between the external vector and S. at an instant of time. 
D = declination of the external vector 

X 

RA = right ascension of the external vector 

x 


Of course a great deal can be known about the behavior of a particular sen- 
sor with respect to an external vector without the complete solution. For instance 
if the momentum vector is known, the minimum angle a sensor ever makes with 
an external vector can be found. From Figure 14 the minimum angle S. makes 
with an external vector is given by 


fl. min = 0, - ARC. 

i ' L i 


if ( I fl L - ARC. 1-61 > o) 


(90) 


a min =0 if (i n L - ARC. | -8 < 0) 


(91) 


where is the angle between the momentum vector and the external vector, 
and 6 is the half angle of any precession which may be present. 
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The time at which 0. min and fi. max occur can also be found rather easily 
in the case where 6=0. From Figure 15 find 


sin D - cos H, sin D, 

X x L L 

= 

sin fl L cos D l 
cos D sin (RA - RA, ) 

» \r X X L 

sin X - : 

sin 


(94) 


(95) 


The time at which a sensor makes a minimum angle with the external vector 
is then 


c a 


= t + 


(0 -X) - AZSEN. 

( 4 > + 


(96) 


where t s is the time at which the aspect sensor crosses the sun, and Q is the 
solution to equations (79) and (80). 

The time at which a sensor makes a maximum angle with the external vector 
is 


tfl. =tfc. . +-r^-r- (97) 

1 max 1 mln (<£ + 0) 


li 6 / 0, the general method of solving for Q. must be used a number of 
times in a roll period and the minimum value of fi. set equal to H. . . 

i i min 

Similarly for Cl. 

^ i max 
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APPENDIX A 
Plots of At. and h 

A 

Solutions of Equation (61) 
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APPENDIX B 

Plots of At A and At B 
Solutions of Equation (63) 
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APPENDIX C 

Plots of At c and /3 
Solution of Equation (68) 
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APPENDIX D 


Tables of Error in Momentum 
Vector Position Caused by 1 ° Error in 
8 for Various Values of 77 and J3. 
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ETA- 10.0 
SETA 


DELTA 

10 

20 

30 

60 

50 

60 

70 

80 

90 

100 

110 

120 

130 

160 

150 

160 

170 

2.0 

1.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

6.0 

1.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

6.0 

1.1 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

8.0 

1.1 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

10.0 

1.2 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

12.0 

1.3 

2.6 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

16.0 

1.6 

2.1 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

16.0 

1.8 

2.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

18.0 

2.6 

2.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

20.0 

6.2 

2.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

22.0 

0 . 

2.2 

3.8 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

26.0 

0 . 

2.5 

3.2 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

26.0 

0 . 

3.0 

2.9 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

28.0 

0 . 

6.5 

2.9 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

30.0 

0 . 

39.6 

2.9 

78.1 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

32.0 

0 . 

0 . 

3.1 

5.1 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

36.0 

0 . 

0 . 

3.5 

6.2 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

36.0 

0 . 

0 . 

6.2 

3.8 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

38.0 

0 . 

0 . 

6.2 

3.7 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

60.0 

0 . 

0 . 

59.5 

3.7 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

62.0 

0 . 

0 . 

0 . 

3.9 

6 . 3 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

66.0 

0 . 

0 . 

0 . 

6 . 3 

5 . 1 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

66.0 

0 . 

0 . 

0 . 

5.2 

6.6 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

68.0 

0 . 

0 . 

0 . 

7.7 

6.6 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

60.0 

0 . 

0 . 

0 . 

79.6 

6.6 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

62.0 

0 . 

0 . 

0 . 

0 . 

6.6 

7.2 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

56.0 

0 . 

0 . 

0 . 

0 . 

5 . 1 

5.8 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

56.0 

0 . 

0 . 

0 . 

0 . 

6 . 1 

5.2 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

58.0 

0 . 

0 . 

0 . 

0 . 

8 . 9 

5.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

60.0 

0 . 

0 . 

0 . 

0 . 

99.2 

5.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

62.0 

0 . 

0 . 

0 . 

0 . 

0 . 

5.2 

8.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

66.0 

0 . 

0 . 

0 . 

0 . 

0 . 

5.7 

6.6 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

66.0 

0 . 

0 . 

0 . 

0 . 

0 . 

6.8 

5.7 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

68.0 

0 . 

0 . 

0 . 

0 . 

0 . 

9.9 

5.5 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

70.0 

0 . 

0 . 

0 . 

0 . 

0 . 
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0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

154.0 

0 . 

0 . 

l . 3 

1.0 

1.0 

1.2 

1.8 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

156.0 

0 . 

0 . 

l . 5 

1.0 

1.0 

1.2 

2 . 1 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

158.0 

0 . 

0 . 

2. 2 

1.0 

1.0 

1.2 

3.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

160.0 

0 . 

0 . 

60.0 

1. 1 

1.0 

l . 2 

140.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

162.0 

0 . 

0 . 

0 . 

1 . 1 

1.0 

1.3 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

u. 

0 . 

0 . 

0 . 

0 . 

164.0 

0 . 

U . 

0 . 

1.2 

1 .0 

1.4 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

166.0 

0 . 

0 . 

0. 

1 .4 

l .0 

1.6 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

168.0 

0 . 

0 . 

0 . 

1.9 

1.0 

2.2 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0. 

0 . 

0 . 

0 . 

0 . 

l 70.0 

0 . 

0 . 

0 . 

8 C .0 

1.0 

1 20.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0. 

0 . 

0 . 

0 . 

0 . 

0 . 

172.0 

0 . 

0 . 

0 . 

0 . 

1.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

l 74.0 

0 . 

0 . 

0 . 

0 . 

l . 0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

176.0 

0 . 

0 . 

0 . 

0 . 

l .0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

l 78.0 

0 . 

0 . 

0 . 

0. 

1.0 

0 . 

0 . 

0 . 

0 . 

0. 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 
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ETA > 1 <, 0.0 

beta 


delta 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

2.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

1.0 

0 . 

0 . 

0 . 

4,0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

1.0 

0 . 

0 . 

0 . 

6.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

1.0 

0 . 

0 . 

0 . 

6.0 

0 . 

0 . 

◦ . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

l.o 

0 . 

0 . 

0 . 

10. 0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

99.8 

1.0 

59.8 

0 . 

0 . 

12.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

1.8 

1.0 

1.5 

0 . 

0 . 

14.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

1.5 

1.0 

1.2 

0 . 

0 . 

16.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

1.4 

1.0 

l.l 

0 . 

0 . 

18.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

1.3 

1.0 

l.l 

0 . 

0 . 

20.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

119.5 

1.3 

1.0 

1.0 

39.7 

0 . 

22.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.5 

1.2 

1.0 

1.0 

1.6 

0 . 

24.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.0 

1.2 

1.0 

1.0 

1.3 

0 . 

26.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

1.8 

1.2 

1.0 

1.0 

1 . 1 

0 . 

28.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

1.6 

1.2 

1.0 

1.0 

1 . 1 

0 . 

30.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

9.2 

1.5 

1.2 

1.1 

1.0 

1.0 

19.6 

32.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

3 . 1 

1.5 

1.2 

1.1 

1.0 

1.0 

1.4 

34.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.4 

1.4 

1.2 

1.1 

1.0 

1.0 

1. 1 

36.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2 . 1 

1.4 

1.2 

1.1 

1.0 

1.0 

1.0 

38.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

1.9 

1.4 

1.2 

1.1 

1.0 

1.0 

1.0 

40.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

10.7 

1.8 

1.4 

1.2 

l.l 

1.1 

1.0 

1.0 

42.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

3.5 

1.7 

1.4 

1.2 

l.l 

l.l 

1.1 

l.l 

44.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.7 

1.6 

1.4 

1.2 

1.1 

1.1 

1.1 

1.2 

46.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.4 

1.6 

1.3 

1.2 

1.2 

1.1 

1.1 

1.4 

48.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.1 

1.6 

1.3 

1.2 

1.2 

1.2 

1.2 

2.1 

50.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.0 

1.5 

1.3 

1.3 

1.2 

1.2 

1.3 

19.9 

52.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

3.9 

1.9 

1.5 

1.4 

1.3 

1.2 

1.3 

1.4 

0 . 

54.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

3.0 

1.8 

1.5 

1.4 

1.3 

1.3 

1.3 

1.6 

0 . 

56.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.5 

1.7 

1.5 

1.4 

1.3 

1.3 

1.4 

2.0 

0 . 

58.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.3 

1.7 

1.5 

1.4 

1.3 

1.3 

1.5 

3.1 

0 . 

60.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.1 

1.6 

1.5 

1.4 

1.4 

1.4 

1.6 

39.9 

0 . 

62.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

4.1 

2.0 

1.6 

1.5 

1.4 

1.4 

1.4 

1.8 

0 . 

0 . 

64.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

3.1 

1.9 

1.6 

1.5 

1.4 

1.4 

1.5 

2.0 

0 . 

0 . 

66.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.6 

1.8 

1.6 

1.5 

1.4 

1.5 

1.6 

2.5 

0 . 

0 . 

68.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.4 

1.8 

1.6 

1.5 

1.4 

1.5 

1.7 

3.8 

0 . 

0 . 

70.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.2 

1.7 

1.5 

1.5 

1.5 

1.5 

1.9 

59.9 

0 . 

0 . 

72.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

4.1 

2.0 

1.7 

1.5 

1.5 

1.5 

1.6 

2 . 1 

0 . 

0 . 

0 . 

74.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

3.1 

1.9 

1.6 

1.5 

1.5 

1.5 

1.7 

2.4 

0 . 

0 . 

0 . 

76.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.6 

1.9 

1.6 

1.5 

1.5 

1.6 

1 .8 

2.9 

0 . 

0 . 

0 . 

78.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.4 

1.8 

1.6 

1.5 

1.5 

1.6 

1.9 

4.5 

0 . 

0 . 

0 . 

80.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2.2 

1.7 

1.6 

1.5 

1.6 

1.7 

2.1 

79.9 

0 . 

0 . 

0 . 

82.0 

0 . 

0 . 

0 . 

0 . 

0 . 

4.0 

2.0 

1.7 

1.6 

1.5 

1.6 

1.7 

2.3 

0 . 

0 . 

0 . 

0 . 

84.0 

0 . 

0 . 

0 . 

0 . 

0 . 

3 .C 

1.9 

1.7 

1.6 

1.6 

1.6 

1.8 

2 .6 

0 . 

0 . 

0 . 

0 . 

86.0 

0 . 

0 . 

0 . 

0 . 

0 . 

2.6 

1.8 

1.6 

1.6 

1.6 

1.6 

1.9 

3.2 

0 . 

0 . 

0 . 

0 . 

88.0 

0 . 

0 . 

0 . 

0 . 

0 . 

2.3 

1.8 

1.6 

1.6 

1.6 

1.7 

2.0 

4.9 

0 . 

0 . 

0 . 

0 . 

90.0 

0 . 

0 . 

0 . 

0 . 

0 . 

2.1 

1.7 

1.6 

1.6 

1.6 

1.7 

2.2 

99.9 

0 . 

0 . 

0 . 

0 . 

92.0 

0 . 

0 . 

0 . 

0 . 

3.8 

2.0 

1.7 

1.6 

1.6 

1.6 

1.8 

2.4 

0 . 

0 . 

0 . 

0 . 

0 . 

94.0 

0 . 

0 . 

0 . 

c . 

2.9 

1.9 

1.6 

1.6 

1.6 

1.6 

1.9 

2.8 

0 . 

0 . 

0 . 

0 . 

0 . 

96.0 

0 . 

0 . 

0 . 

0 . 

2.4 

1.8 

1.6 

1.5 

1.6 

1.7 

2.0 

3.4 

0 . 

0 . 

0 . 

0 . 

0 . 

98.0 

0 . 

0 . 

0 . 

0 . 

2.2 

1.7 

1.6 

1.5 

1.6 

1.7 

2.1 

5.2 

0 . 

0 . 

0 . 

0 . 

0 . 

100. 0 

0 . 

0 . 

0 . 

0 . 

2.0 

1.6 

1.5 

1.5 

1.6 

l.H 

2.3 

119.9 

0 . 

0 . 

0 . 

0 . 

0 . 

102.0 

0 . 

0 . 

0 . 

3.4 

1.8 

1.6 

1.5 

1.5 

1.6 

1.8 

2.5 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

104.0 

0 . 

0 . 

0 . 

2.6 

1.7 

1.6 

1.5 

1.5 

1.6 

1.9 

2.8 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

106.0 

0 . 

0 . 

0 . 

2.2 

1.6 

1.5 

1.5 

1.5 

1.7 

2.0 

3.5 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

108.0 

0 . 

0 . 

0 . 

1.9 

1.6 

1.5 

1.5 

1.5 

1.7 

2.1 

5.3 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

110.0 

0 . 

0 . 

0 . 

1.8 

1.5 

1.5 

1.5 

1.6 

1.7 

2.3 

139.9 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

112.0 

0 . 

0 . 

3.0 

1.7 

1.5 

1.4 

1.5 

1.6 

1.8 

2.5 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

114.0 

0 . 

0 . 

2.2 

1.6 

l .4 

1.4 

1.5 

1.6 

1.9 

2.8 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

116.0 

0 . 

0 . 

1.9 

1.5 

1.4 

1.4 

1.5 

1.6 

1.9 

3.5 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

118.0 

0 . 

0 . 

1.7 

1.4 

1.4 

1.4 

1.5 

1.6 

2.0 

5.2 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

120.0 

0 . 

0 . 

1.5 

1.4 

1.3 

1.4 

1.5 

1.7 

2.2 

159.8 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

122.0 

0 . 

2.4 

1.4 

1.3 

1 . 3 

1.4 

1.5 

1.7 

2.4 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

124.0 

0 . 

1.8 

1 . 3 

1.3 

1 . 3 

1.4 

1.5 

1.8 

2.7 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

126.0 

0 . 

1.5 

1.3 

1.3 

1 . 3 

1.4 

1.5 

1.8 

3.3 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

128.0 

0 . 

1.4 

1.2 

1.2 

1.3 

1.3 

1.5 

1.9 

5.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

130.0 

0 . 

1.2 

1.2 

1.2 

1.2 

1.3 

1.5 

2.0 

177.4 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

132.0 

1.6 

1.2 

l.l 

1.2 

1.2 

1.3 

1.6 

2.2 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

134.0 

1.3 

1.1 

1.1 

1.2 

1.2 

1.4 

1.6 

2.5 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

136.0 

1 . 1 

1. 1 

1.1 

l.l 

1.2 

1.4 

1.7 

3.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

138.0 

1.0 

l.o 

l.l 

1.1 

1.2 

1.4 

1.7 

4.5 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

140.0 

1.0 

1.0 

1.1 

l.l 

1.2 

1.4 

1.8 

159.9 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

142.0 

1.0 

1.0 

1.0 

l.l 

1.2 

1.4 

2.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

144.0 

1.1 

l.o 

1.0 

l.l 

1.2 

1.4 

2.2 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

146.0 

1.2 

1.0 

1.0 

l.l 

1.2 

1.5 

2.7 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

148.0 

1.7 

1.0 

1.0 

1.1 

1.2 

1.5 

3.9 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

150.0 

19.9 

L.L 

1.0 

1.1 

1.2 

1.6 

140.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

152.0 

0 . 

1.1 

1.0 

1.0 

1.2 

1.7 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

154.0 

0 . 

1.2 

1.0 

1.0 

1.2 

1.9 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

156.0 

0 . 

1.4 

1.0 

1.0 

1.2 

2.2 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

158.0 

0 . 

2.0 

1.0 

1.0 

1.2 

3.2 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

160.0 

0 . 

40.0 

1.0 

1.0 

1.3 

120.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

162.0 

0 . 

0 . 

l.l 

1.0 

1.3 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

164.0 

0 . 

0 . 

1.2 

1.0 

1.4 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

166.0 

0 . 

0 . 

1.3 

1.0 

1.6 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

168.0 

0 . 

0 . 

1.8 

1.0 

2. 2 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

170.0 

0 . 

0 . 

60.0 

1.0 

100. 0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

172.0 

0 . 

0 . 

0 . 

1.0 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 
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APPENDIX E 

Removing Uncertainty in t 0 Caused by a 
Quantized Time Base in the Spacecraft 

Assume that 8 ~- r - 0. 

Assume that the time base is quantized in the spacecraft such that the time 
t 0 at which the digitized fan crosses the sun is known to lie in the time interval 
between 1 0B and t QB . One revolution later the next time the fan crosses the 
sun the "see sun" time is tj which lies in the interval between t 1A and t JB . 
See Figure A. 
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Figure A— Sun sightings defined within a quantized time base 
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We can also write 




t . > (t 1A - T) t 1 (t 1B - T) 

min v 1A max v IB 

t > (t, A - 2T) t £ (t, R - 2T) 

min v 2A max 2B 

t £ <t,. -3 T) t £ (t, R - 3T) 

min x 3A max 3B 

. • 


It should be possible to find a t and a t which will satisfy all the 

r min max J 

above inequalities. If the inequalities are satisfied over a sufficiently large 
number of revolutions, the value of t . should approach that of t . Thus it 
is possible to define the value of t 0 within any desired accuracy by merely 
solving the inequalities over a sufficiently large number of revoltuions. 

This method of solution is vulnerable to noise. It is therefore advisable to 
test to make sure that t . £ t holds at all times. 

min max 
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APPENDIX F 


Horizon sensors which operate in the visible spectrum may be used as long 
as the earth as seen from the satellite is fully illuminated. From Figure AA it 
can be seen that the spacecraft sees a fully illuminated earth as long as 


V + (90° - p) < 90° 


or 


V < P 



Figure AA-Relationships of spacecraft to sunlit earth 
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